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Abstract 

cp ; 

In the Optimal Velocity Model proposed as a new version of Car Following Model, 
it has been found that a congested flow is generated spontaneously from a homoge- 
neous flow for a certain range of the traffic density. A well-established congested flow 
obtained in a numerical simulation shows a remarkable repetitive property such that 
^ ■ the velocity of a vehicle evolves exactly in the same way as that of its preceding one 

except a time delay T. This leads to a global pattern formation in time development 
of vehicles' motion, and gives rise to a closed trajectory on Ax-v (headway- velocity) 
plane connecting congested and free flow points. To obtain the closed trajectory an- 
alytically, we propose a new approach to the pattern formation, which makes it pos- 
sible to reduce the coupled car following equations to a single difference-differential 
equation (Rondo equation). To demonstrate our approach, we employ a class of 
linear models which are exactly solvable. We also introduce the concept of "asymp- 
totic trajectory" to determine T and vb (the backward velocity of the pattern), the 
global parameters associated with vehicles' collective motion in a congested flow, in 
terms of parameters such as the sensitivity a, which appeared in the original coupled 
equations. 
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1 Introduction 



Traffic flow is one of the most interesting phenomena of many-body systems which may 
be controlled by basic dynamical equation. Recent developments in study of traffic flow 
has brought a renewed interest in microscopic approaches, such as Optimal Velocity Model 
(OV-model) |], @, [J, which is a new version of Car Following Model ||, |^, [J, Cellular 
Automaton Models 0, || || , Coupled Map Lattice Models JTOf and Fluid Dynamical Models 



nj. The OV-model, among others, has especially attracted interest because it provides 
us with a possibility of unified understanding of both free and congested traffic flows from 
common basic dynamical equations. Unlike traditional Car Following models, it introduces 
Optimal Velocity Function V(Ax) as a desirable velocity depending on headway distance 
Ax. The basic equation of OV-model for a series of vehicles on a circuit of length L is, 

x n {t) = a{V{Ax n (t)) - x n {t)} n = l,2---N, (1) 

where x n denotes the position of the n-th vehicle, Ax n = x n -\ — x n headway and N the 
total vehicle number. The constant parameter a is the sensitivity. A driver accelerates (or 
decelerates) his vehicle in proportion to the difference between his velocity and the optimal 
velocity V(Ax). As easily noticed, a homogeneous flow is a solution to Eq. (|l|). In such 
a flow, vehicles have a common headway L/N, which is the inverse of the vehicle density. 
Stability of homogeneous flows is analyzed within a linear approximation jl], 0; it is stable 
for / = V'(L/N) < f c and unstable for / > f c . The critical value is found to be f c = a/2. 

In order to demonstrate that the OV-model describes "spontaneous generation of con- 
gestion", numerical simulations were made using Eq. ([!]). It was found that for / < f c , i.e., 
if the density is above the critical value, a slightly perturbed homogeneous flow develops to 
a congested flow after enough time. The congested flow consists of alternating two distinct 
regions; congested regions (high density), and smoothly moving regions or free regions (low 
density). In this way the traffic congestion occurs spontaneously in the OV-model. This 
phenomenon can be understood as a sort of phase transition from a homogeneous flow 
state to a congested flow state [|I], 0. 

A remarkable feature of the well-established congested flow is that the velocity of the 
n-th vehicle x n has the same time dependence as that of the preceding ((n — l)-th) vehicle 
except a certain time delay T. It is also found that the global pattern moves backward 
with a velocity vb- This kind of behavior of vehicles may be called "repetitive pattern 
formation". It leads to formation of a closed trajectory ("limit cycle") on Ax-v plane, 
along which representative points for all the vehicles move one after another. The con- 
vergence of vehicles' trajectories to a closed trajectory signals the congestion in a traffic 
flow. Therefore the determination of the closed trajectory is one of the most important 
subjects to understand congested flows. However this has been done mostly in computer 
simulations. The purpose of the present paper is to obtain this closed trajectory directly 
by an analytical method. Actually it has been found that, in the vicinity of the critical 
point for the congested flow, equation ([!]) can be reduced to the modified KdV equation by 



the dynamical reduction method [|12|]. In this paper we propose a new analytical approach 
to the pattern formation, which may be applicable to any congested flow. 

We argue that, once the repetitive pattern is formed, the coupled car following equations 
reduce to a single difference-differential equation (Rondo equation) for a universal function 
(Rondo function). A Rondo function determines a closed trajectory on Ax-v plane. To 
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make the Rondo equation tractable, we have simplified our question on the following two 
points: firstly we have assumed that OV-functions are piece-wise linear; secondly we have 
concentrated our attention to an asymptotic trajectory, which is the key concept to be 
explained in the next section. We would like to stress that our method does not lose 
its generality by making the above assumption on OV-functions: an OV-function to be 
obtained from real data may be approximated by a piece-wise linear function. 

With the above simplifications, we have solved the Rondo equation for each model 
and given an asymptotic trajectory on the Ax-v plane. Our result clearly tells us that, 
once an OV-function and the sensitivity a are given, an asymptotic trajectory is uniquely 
determined; this then implies that the parameters T and vb for a collective motion of 
vehicles are given as a function of a. Therefore our approach provides us with a method 
to determine a dependence of the global parameters T and vb- 

This paper is organized as follows. Next section summarizes the main results obtained 
from numerical simulations of the OV-model, with emphasis on the pattern formation in 
a congested flow. The concept of an asymptotic trajectory is explained in this section. As 
will become clear in later sections, an asymptotic trajectory is a very important concept 
to understand the OV-model. We derive the difference-differential equation and present a 
general strategy on how to solve it in section |^. This part summarizes the central idea of 
this paper. In order to demonstrate how the Rondo approach works, we analytically solve, 
in section f|, some simple models with piece-wise linear OV-functions. Our first model has 



been investigated by Sugiyama and Yamada[13|. Here we solve this model in the context of 



the Rondo approach. Although an asymptotic trajectory is very close to a real trajectory 
observed in a computer simulation, it is not exactly the same to the latter. We describe 
some aspects of real trajectories based on our knowledge of the asymptotic one in section |^. 
The final section is devoted to summary and discussions. 



2 Pattern Formation in OV-Model 

Let us recollect what we have learned with numerical simulations of an OV-model |], |3| ■ 
Suppose a simulation is performed with a given OV-function and a fixed sensitivity a. 
After a congested flow is well-established, typical features of the repetitive behavior can 
be observed in the following two figures. 

Figs. |l](a) and (b) show that vehicles move in alternating regions of free and congested 
flows. It is recognized that every vehicle behaves in the same manner as its preceding one 
with a certain time delay T: as a result, the congested region moves backward with the 
velocity vb- Once the location of any vehicle, say, the n-th vehicle, is given as a function 
of t, we may reproduce the pattern in Fig. |l|(a) by plotting a series of functions shifted in 
time and position by T and vbT appropriately. Therefore we expect that a congested flow 
may be completely determined by a function of t and global parameters T and VbT. The 
precise specifications of our approach to this repetitive behavior will be explained in the 
next section. 

Fig. |l|(b) clearly shows there exists a "limit cycle" on the Ax-v plane, a closed curve 
with two cusps at points C (Axc,vc) and F (Axf,vf), both of which are on the OV- 
function. At these cusp points, we find V'(Ax) < a/2, which means that homogeneous 
flows with such headways and velocities are linearly stable. Representative points for all 
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Time Headway 

Figure 1: A result of the numerical simulation with 20 vehicles on a circuit with the 
circumference L = 40. The sensitivity and the OV- function are chosen as a = 1.0 and 
V(h) = tanh(/i — 2) + tanh(2). (a) Trajectories of vehicles passing through a congested 
region, (b) The "limit cycle" on the Ax-v (headway- velocity) plane. The cusp points C and 
F correspond to congested and free regions respectively. Representative points for vehicles 
move anti-clockwise along this "limit cycle" . They run fast on curves connecting the two 
cusps, while they stay around cusps for a while. The thin line shows the OV-function. 



vehicles move on this "limit cycle" in the same direction as anti-clockwise. It follows from 
the conservation of flow, 

Ax c ~ Ax F ~ T' U 
that the straight line connecting C and F has the slope T _1 and intersects with the vertical 
axis at — i>b-[] 

In the rest of this section, we would like to explain the concept of an asymptotic 
trajectory on the Ax-v plane. Suppose that a vehicle passes through two congested regions 
different in size on a circuit. Then the representative point moves on curves as shown 
in Fig. |2|. Each trajectory may not form an actual cusp, rather it will form a round 
shaped tip. Also we find that the larger the congested region, the sharper the shape 
of the tip: the minimum velocity of the vehicle is smaller for a longer congestion. It is 
rather easy to imagine that for a very short congestion the vehicle cannot decelerate itself 
enough to reach the velocity appropriate for a longer congested region. If we plot minimum 
velocities for longer and longer congested regions, we would find a limiting value for the 
minimum velocities. This value must be realized for an infinitely long congested region. 
With a similar argument we find the limiting value for maximum velocities corresponding 
to an infinitely long free region. We may imagine the following extreme situation: a 
vehicle starting from an infinitely long free (congested) region goes toward an infinitely 
long congested (free) region. The trajectory for this limiting situation will be called as 
a decelerating (an accelerating) asymptotic trajectory. Combining them we would find a 
closed curve with two real cusps on the OV-function. 

^his relation is an approximate one except for an asymptotic trajectory to be discussed below. 
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Figure 2: Trajectories for a vehicle passing through two congested regions different in 
size. The outer curve corresponds to the larger congested region. The thin line is the 
OV- function. Parameters are the same as those in Figure 1. 

The duration for a vehicle to stay in a congested region would obviously get longer for a 
larger congested region. For an asymptotic trajectory it becomes infinite. This does fit to 
our linear analysis since the behavior of a vehicle is controlled by exponential functions in 
time. In section |5|, we shall see that exponential functions determine a curve near a cusp. 

3 Rondo Approach 

We begin our description of the Rondo approach with two basic assumptions: 

1. Velocities of the n-th and (n — 1) -th vehicles have exactly the same time dependence 
if a certain time delay T is taken into account, x n _i(t) = x n (t + T). 

2. The pattern of traffic flow moves backward with a constant velocity v B - 
The above described properties are expressed as, 

x n -i{t) = x n (t + T) + v B T. (3) 

All the vehicles' behavior is represented with a single universal function F(t) = £.„(?;) ;0 

x n _ fc (i) = F{t + kT) + kv B T. (4) 

With n-th vehicle's headway given as 

Ax n (t) = F(t + T) - F{t) + v B T, (5) 

N coupled car following equations (|l|) are reduced to a single difference-differential equation 
for F(t), 

-F(t) + F(t) = V(F(t + T) — F(t) + v B T). (6) 
a 

2 This assumption bares some similarities with the travelling wave ansatz (<fi(x,t) — > f(x — vt)) for the 
wave equation. 
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In the following it will be called the Rondo equation. 

In this paper, we will seek the Rondo function F(t) for the asymptotic trajectory.^ 
Before studying concrete models, let us consider its generic properties. Since the posi- 
tion and the velocity of vehicles are obviously continuous in time, F(t) is a continuously 
differentiable function. 

An asymptotic trajectory connects the points F and C, each of which corresponds to an 
infinitely long free or congested region (an approximately homogeneous flow) satisfying the 
stability condition mentioned in section |^. Therefore F(t) should be homogeneous flows 
asymptotically in the infinite past and future: F(t) — > const, as t — > ±00. An asymptotic 
trajectory interpolates two stable solutions of the equation ([!]). In this sense F(t) may be 
regarded as a "kink solution". 

Like OV-models studied in earlier papers, each model in the next section has an OV- 
function which is symmetric with respect to a point, S(Aa;s, t>g).Q So we assume this 
property in the following arguments and quote our result in appendix [A]. Two end points 
of an asymptotic trajectory, C^Axc^c) an d F(Ax F ,v F ), are symmetric with respect to S. 
Three points C, F and S are on a straight line with a slope T~ l and an intersection — v B - 
Note that once the slope is given, the intersection is uniquely determined since the point 
S must be on the line. 

As shown in appendix [A], accelerating and decelerating asymptotic trajectories are 
symmetric with respect to S. Therefore it is sufficient to study one of them; in the rest 
of this paper, we will take an decelerating asymptotic trajectory. Here we summarize 
conditions which should be satisfied by the function F(t): 

1) F(t) and F(t) are continuous for any t (F(t) £ C 1 ); 

2) v F = lim F(t) and v c = lim F(t); 

3) v F + v c = 2v s and Ax F + Ax c = 2Ax s , where Ax F = lim (F(t + T) - Fit) + v B T) 

t~> — OC 

and Ax c = lim (F(t + T) - F(t) + v B T). 

t — >+oo 

We would like to explain a way to solve the Rondo equation, which contains an OV- 
function and a sensitivity a as well as T and v B associated the pattern formation. (1) 
First we give the parameter T. By drawing a straight line through the point S with the 
slope T -1 , we find the intersection — v B and the points C and F. (2) Now a is the only 
free parameter of the Rondo equation. If we could solve the equation, we would obtain 
a one-parameter family of (or a-dependent) Rondo functions. (3) Among this family, the 
right Rondo function is selected by requiring that it connects the points C and F. This 
condition also determines a unique value for a. Accordingly we find the a-dependence of 
T. 

4 Piece-wise Linear Function Models 

We consider here a class of models with piece-wise linear OV-functions. The Rondo equa- 
tion is now linearized for all regions of Ax, and therefore exactly solvable. 

3 We will discuss more realistic situations with finite congested regions in section 5. 
4 The symmetry of an OV-function is absolutely not necessary to solve a system in the Rondo approach. 
In the last section, we discuss how to solve the Rondo equation for a generic situation. 
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4.1 Step Function Model 

The first model has the step function for the OV-function, 



V(Ax) = { ° Ax < Axs ( regi ° n ^ (7) 

\ V Ax > Ax s (region II). K } 



This OV-model has been solved in Ref. []n|. Here we explain how this model can be solved 
in our Rondo approach. In this model, the Rondo equation is given by 

-F(t) + F(t) = V 9 (F(t + T) - F{t) + v B T - Ax s ) , (8) 
a 

where 6(x) is the Heviside function. 

In the motion corresponding to a decelerating asymptotic trajectory, the representative 
point for a vehicle moves from region II into region I. Let us take the time t such that the 
point moves into the region I at t — 0, which implies Ax n (0) = Ax$. 

The equation of motion is 

'fW + fwJI" («<?) (9) 



a 

The general solutions for two regions are 



(t > 0). 



F{t) - { C 2 e~ at (t > 0). (10) 

The integration constants are determined as C\ = and C% = Vo from requirements that 
F(t) be a continuous function and asymptotically constant for t — > ±oo. The continuous 
function F(t) is 

r v t (t < o) 

F ® = \ ^(l-e-) (t>0). 

Here we choose the origin of position coordinate so that F(0) = 0. 

The function F(t) with its relation to the headway, Ax n (t) = F(t + T) — F(t) + vbT, 
completely determines a decelerating asymptotic trajectory on the Ax-v plane. From the 
condition 3) in the section |3], the asymptotic trajectory connects symmetric points on the 
OV-function. This gives us a condition Ax n (+oo) + Ax n (— oo) = 2Ax$ = 2Ax n (0), 



vbT + (V + v B )T = 2 



\l-e~ aT ) + v B T 
L a 



(12) 



This may be expressed as the transcendental equation for p = aT 

1 



P + -p-l = 0. (13) 



Since the p is found to be the constant (1.59362..), we obtain 

aT = p = 1.59362.. . (14) 
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Figure 3: (a) The OV-function (thick line) and a decelerating asymptotic trajectory (thin 
line with arrows) for the Step Function Model, (b) The position functions x n (t) and x n -i(t) 
for the decelerating asymptotic trajectory. 



This gives us the a-dependence of T, which was first obtained in Ref. [T3[ . 

It is instructive to see a relation between the function F(t) and the decelerating tra- 
jectory depicted in Fig. |3|. Two curves in Fig. |3](b) correspond to the (n — l)-th and n-th 
vehicles' locations. At t = the n-th vehicle's representative point moves into the region I 
and F(t) is described by an exponential function. Before that time, the function F(t) is 
linear in t. The curve for the (n — l)-th vehicle changes from the linear to the exponential 
behavior at t = —T. It is given via a parallel displacement by the vector (—T,vbT) from 
the curve F(t). For the time t < — T, curves are two parallel straight lines, which implies 
the headway of the n-th vehicle does not change till that time from the infinite past. The 
n-th vehicle has the constant velocity Vo for t < 0. This tells us that the n-th vehicle is in 
the free region for t < — T, indicated by the point F in Fig. 0(a). It is also easy to observe 
that at t = —T the (n — l)-th vehicle starts to decelerate. As a result the headway of the 
n-th vehicle decreases; at t = it reaches the value Axs and the n-th vehicle starts to 
decelerate itself. 

In this model, the points F and C are both characterized as points which are reached 
in the infinite future or past. As a traffic flow, we are describing a soliton like solution 
connecting a half infinite vehicles, running with the velocity Vp and the headway Axf, and 
another half infinite vehicles going into the congested region associated with the point C. 



4.2 Single Slope Function Model 

In the step function model, the Ax-dependence of the Rondo equation is too simple; In- 
dependence is not explicit. So we would like to consider a slightly improved model, whose 
OV-function shown in Fig. f|has a finite slope. Thus we call this model as the Single Slope 
Function Model. The OV-function is characterized by the following parameters: /, the 
slope; Vo, the maximum velocity; Axs, the headway for optimal velocity Vq/2. From the 
linear analysis in the paper fl| , we know that a homogeneous flow becomes unstable and a 
congested flow is expected for 2f > a. This condition is assumed in our analysis here. 



8 




Figure 4: The OV-function for the Single Slope Function Model. / is the gradient in the 
region II. The function is symmetric with respect to the point (Axs, Vq/2). 



The OV-function has sharp bends at 

Ax A = Ax s - \ (15) 

Ax B = Ax s + Yf ( 16 ) 

which divide Ax into three regions, I, II and III, as indicated in Fig. |j. 

We would like to find a Rondo function F(t) for a decelerating asymptotic trajectory, 
along which headway of a vehicle monotonically decreases from Axf to Axo We assume 
it reaches Axb at t — —t and Axa at t — 0. The Rondo equation takes of the form, 



1 



a 




(t>0) 



-F(t) + F(t) = f (F(t + T) — F(t) -5) (-r < t < 0) , (17) 



(t < -r) 



subject to the conditions, 



F(T)-F(0)+v B T = Ax A , (18) 
F(-t + T) -F(-t) +v b T = Ax B . (19) 

Here 5 = Axa — v B T. Note that the time — r is not a free parameter. Rather, it should 



be determined from (|T^) via solving the Rondo equation. 

For regions I and III, the Rondo equation becomes the same as that in the step function 
model. Therefore, we obtain 

p. v _ J const, x e~ at (for region I) , , 

{ ) ~ \ V (for region III). [ ' 

Our purpose in this subsection is to describe a method to find the Rondo function in the 
region II which correctly interpolates those in (EDf). To this end, we may introduce a series 
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of functions for each time interval: 



F(t) 



f F\t) (0 < t) 

Fi(t) (-T<t<0) 

F 2 (t) (-2T <t< —T) 

F m (t) (t < -r) 



(21) 



The first function F l (t) = Fo(t) and the last one F m (t) are for the region I and III, 
respectively. 

Let us find F l (t) ( for t > ). It follows from (EUp that 



F\t) = u e~ at , 
F \ t ) = ^ (l - e~ at 

where we fix again -F^O) = 0. The n-th vehicle's headway is then given by 

Ax n (t) = F\t + T) + v B T - F\t) = ^(1 - e- aT )e- a * + Wb T. 

a 



Since the condition (|18D determines the constant u 

Ax a — vbT a5 

U ° = a l- e -aT 

we find the Rondo function for t > to be 

6 



1-e 



-aT' 



F\t) 



1-e 



-at 



1 - e~ aT 

This leads to a linear relation between the velocity and the headway, 



1-e 



-aT 



(Ax n - v B T). 



(22) 
(23) 



(24) 



(25) 



(26) 



(27) 



In the t — > oo limit, we find that i n — >• and Ax n — > VbT. 

Now we consider Fi(t) ( for —T < t < ). In this time interval, the Rondo equation 
for Fi(t) is expressed as 



-Fx(t) + A(t) = /(F (t + T) - - 5). 



(2f 



in terms of the Rondo function for the region I, F J (t) = -Fo(^)- By using the differential 
operator, 

1 rf2 + 14 + 1, (29) 



V 



af dt 2 f dt 



we may rewrite the above equation as 

VF 1 (t) = F (t + T)-5. 



(30) 
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The general solution may be written as a sum of a particular solution and the solution to 
the homogeneous equation T>F^ ora {t) = 0. It is easy to see that F Q (t + T) — 5 is a particular 
solution, since the function F Q (t) = F J (t) given in ( ^6|) satisfies 

VF (t) = F {t). (31) 

The exponents 7 for a homogeneous solution are 



a I a 2 a , 

T = __ ±i y a/ __ = __ ± ^. ( 32) 

A general solution is given as 

Fi(t) = F (t + T) - 5 + e~% f (Asinujt + B cosiut) . (33) 

The constants A and S are determined as A = a5/u> and B = 0, by the requirement that 
-Fo(t) and Fi(t) must be continuous up to the first derivative at t = 0. Therefore 

Fl (t) = r -^ f (l - e-^+ T )) - 5 + ^e-f* sin^t, (34) 

and the headway and the velocity for the n-th vehicle is given as 

Axjt) = F (t + T) -FUt) +v B T = Ax A - — e'^smcut, 

us 

, . • . . aS n , a5 a, ( a \ , 

x n {t) = Fi(t) = — -e H e a* fu; cos cut - -smwtl . (35) 

We would now like to give general formula for Fk(t) ( for t £ I k = [—kT, —{k — 1)T] ) 
with k > 1. Suppose that -Ffc(t) for t e is known to us and we are trying to find i 7 / c+1 (t) 
for t E Ik+i- Fk+i(t) satisfies the second order linear differential equation, 



with boundary conditions; 



VF k+1 {t) = F k (t + T)-S, (36) 



F k+l {-kT) = F k {-kT), 

F k+1 (-kT) = F k (-kT), (37) 



while F k (t) satisfies a similar equation, 

VF k (t) = F k - 1 (t + T)-5. (38) 

The function F k (t) describes the behavior of the n-th vehicle only for t 6 I k . However 
we will find it useful to define the function by the relation (j38|) even outside the interval 
ifc. The function used outside the interval will be denoted as F k (t). The difference of ( |36| ) 
and ( |38|) gives us the following equation for t G I k +i = [— {k + 1)T, — kT], 

v(F k+l {t)-F k {t))=F k {t + T)-F k „ l {t + T) (k>l). (39) 
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From ( p0|) and (J3T[) we find 

P - F (t)) = F (t + T) — F (t) -5 = 5 (e~ at - l) , (40) 

for teh = [-T,0]. 

By using the function G k (t) defined in the relation, 

F k+1 (t) = F k {t) + G k+1 (t + kT), (41 ) 

9|) and (|40|) are rewritten into the following equations for — T < t < 0, 

VG k+1 (t) = G k (t) (k > 1), (42) 
DG 1 (t) = (5( e - at -l) =G {t). (43) 



Note that Go(t) defined in fl43| ) satisfies 2Xro(t) = Go(t). The conditions fl3"T|) become 
G fc (0) = and G k (0) = ( for k > 1 ). 

The Rondo function F(t) for the region II is given as a sum of G k (t), 

F ll (t) = F (t) + £ - A;T)G fc+1 (t + fcT), (44) 

fc>0 

for t > —t. In the appendix [B] we will give general solutions to differential equations fl42"D 
and (P^). from the appendix, 

dft) = G (t) + — e-t*sinu;t = 5fe- at -l) + — e"** srau;*, (45) 
is consistent to (|34]) . Similarly G 2 (i) is given as 

G 2 (t) = Gift) - —— —7— e~§* {cutcoscut - sincut} . (46) 
uj 4/ — a 

So F 2 (t) for tel 2 = [-2T, — T] becomes 

F 2 (t) =F 1 (t) + G 2 (t + T) 

= FAt) + 5 ( e -^ +T ) - l) + -e~f C +T ) sincu(t + T) 



ab 2/ _|( t+T ) 



e -3l*+'J-J + T) cosc;(t + T) - sina;(t + T)} , (47) 



u 4/ - a 

while the headway for the n-th vehicle is 

Ax n (t) = Axa e _ 2*sina;t 

+ ^J±l_ e -^T) ( u u + T ) C0SLc ,( t + T \ _ sinw ( t + r u . ( 48 ) 
4/ — a 

The general formula in appendix || may be used further to generate F 3 , F 4 , • • •, needed to 
describe the trajectory in the region II. 
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Figure 5: Sensitivity dependence of T and r for the Single Slope Function Model. The 
condition, T > r, is assumed, which is satisfied for a < 0.98857.. . The dotted line shows 
the a-T relation for the Step Function Model. 



Let us consider the region III. At t — — r, the function fl44|) for the region II must be 
continuously connected to F (t) including their first derivatives. This condition yields 
F m (t) for t < -r as 



F m (t) = V (t + r) + F n (-r). 
The continuity condition for the first derivative requires also that 



F u ( 



r-III/ 



Vn. 



(49) 
(50) 



This completes our construction of an asymptotic trajectory. 

Now let us find the a-dependence of T for the present model. The time — r may be 
determined by the condition (|19"D , F(— r + T) — F n (— r) + vbT = Axb] then Eq. fl5Cf) gives 
us a relation between a and T. By using ([0]) for F ll (— r), concrete expressions of ( |I9"D and 
( |50D may be obtained. For — T < — r, F n (— r) is simply given by r) and the above 
conditions are expressed as: 



a<5 a . 
— e 2 sin u;r 

UJ 



ad 



aT 



-e aT + 



a5 



Ax B - Ax A = (51) 

-f 7 leu cos ut H — sin cur] = Vq. (52) 
e"- 1 — 1 cu V 2 / 

The requirement of symmetry, Ax n (+oo) + Ax n (— oo) = 2Axs, mentioned in section |] 
becomes vbT + (Vo + vb)T = 2Axs- This relation and Eq. (|16|) allow us to express 
5 = Ax a — VbT as, 

Vo 



6=(JT- 1) 



(53) 



Finally, we reach to the coupled equations which determine T and r for given slope / and 
sensitivity a: 

(fT — l)e2 r sin cur 



2cu 



(e- -l)[(/-f) 



a 

sin cur 



cu cos cur 



(54) 



cue 2 ' 
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where tu is given in Eq. (p2|) . For / = 1.0, Eq. (|54|) is solved numerically to give the 
a-dependence of T and r as shown in Fig. |^. Since we used Fi(t) for F n (t), Eq. ( |54"D is 
valid only for r < T (r coincides with T when a = 0.98857.. at r = T = 1.74027..). 

We observe in Fig. [5] that T behaves as 1/a for small a: this implies that, for congested 
flows to be formed, the delay T must be larger for less sensitive drivers. In the limit of 
/ — > oo, the present model reduces to the Step Function Model, in which aT = p and 
t = 0. This may be confirmed with Eq. (|54"D since aT reaches a constant, p = 1.59362.., 
and ar behaves like 2a/ (fp) when a/ f goes to zero.0 



a=1.8 






/ / / 






/ 

t = -T/ / 




t = o j// 
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Figure 6: Thin curves show "limit cycles" for the Single Slope Function Model obtained 
from simulations with (a) a = 1.8, (b) a = 1.5, (c) a = 1.0 and (d) a = 0.9. Trajectories 
generated from Fi(t) for t G [— T, 0] are drawn with thick curves. The curves inside the 
region II (1 < Ax < 3) are parts of asymptotic trajectories. The dotted line shows the 
OV-function with / = 1.0, V = 2.0 and Ax s = 2.0. 

In order to see the validity of our approach, let us compare our results and trajectories 
obtained by simulations. In Fig. || thick curves show parts of the asymptotic trajectories, 

5 Note that Eq. (M) may be rewritten in terms of rescaled variables aT, ar and a/ f. 
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to be determined by the function F%(t), for a = 1.8, 1.5 (r > T), a = 1.0 (r ~ T) 
and a = 0.9 (r < T). These curves are to be compared with numerical simulations 
shown as thin curves. The function Fi(t) is enough to give an asymptotic trajectory for 
a < 0.98857.., as mentioned above. It is expected that when a gets closer to its critical 
value (a cr iticai = 2/ = 2.0, in this case), we need functions Fk(t) with higher k to form an 
entire trajectory. 

There appears a flat trajectory again in the region III. As in the Step Function Model, 
it takes time T for a vehicle to move on the flat trajectory and there is only one vehicle 
traveling on this interval. 

4.3 Double Slope Function Model 

The OV-function for the single slope model has flat regions I and III like the step function 
model. For those regions the Rondo equation does not depend on F{t + T): a vehicle does 
not react to the motion of the preceding one. This motivates us to consider a more realistic 
model with an OV-function which has non-zero gradient for any headway. 



I 


II 


III 








' fx 




fl 















Axa Ats Axb At 



Figure 7: The OV-function for the Double Slope Function Model, which is symmetric 
with respect to the point (Axs, V(Aa;s))- R has the gradient fi (^2) in the regions I and 
III (II). 

The OV-function has a slope f\ in regions I and III, and f'2 in region II (see Fig. 0): 

{fiAx Ax < Ax a (region I) 

f 2 [Ax - (1 - k)Ax a } Ax a < Ax < Ax B (region II) (55) 

fi [Ax + (kT 1 - l)(Ax B - Ax A )] Ax b < Ax (region III). 

where k = fx/ fz- Obviously this function is symmetric around (Axg, V(Axs)), where Ax$ 
is the middle point of the region II. Here we should note that the sensitivity a must satisfy 
fi < a/2 < f 2 for generation of the congestion in this model, since the homogeneous flow 
is expected to be linearly unstable only in the region II. 

As in the previous subsection, we assume that the headway reaches Axb at t = — r and 
Axa at t = 0. The solution of the Rondo equation (|6|) must satisfy the conditions (pT8|), 
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(19). The Rondo function F(t) for three regions will be denoted as follows; 

f F\t) (0 < t) 

F{t) = I F u {t) {-r<t< 0) (56) 

{ F lll (t) (t < -r). 

For -F : (t) (t > 0), the Rondo equation is given by 

PiF ! (t) = F r (t + T) + Vi jT, (57) 

where 

^ Id 2 Id, 

B -^ + AS + L (58) 

To find a solution to the homogeneous equation T>iF^ om {t) = Fy\ ora (t+T), we use the ansatz 
F^ om (t) ~ e 7 * which gives an equation for the exponent 7, 



~ 2 +7 = /i(e^ T -1). (■-»») 



m 



a 

As long as the condition f\ < 1/T holds|] there are two real solutions: the negative —7 
and the positive one 7 out . Because of the asymptotic behavior, F(t) — > vq as t — > +00, 
only the exponent — 7; n is relevant to the function F l (t). Adding a particular solution of 
(|57D, we obtain the solution subject to conditions ( |T8D and F(0) = as 

F ! (t) = Vc t + (Ax A - Ax c ) }~ 6 _T T - (60) 

1 — e 'in- 1 

By calculating Ax n (t) and x„(t) from (pO), 

Ax n (t) = Ax c + (Ax A - Ax c )e^\ (61) 

£ n (t) = + (Ax A - Axc) - 7i ° T e~ 7 '" f , (62) 

1 — e 7in 

we obtain a linear trajectory given by 

x n (t)-v c = - — 7l \ r (Aar w (t) - Ax c ) . (63) 
1 — e _w 

In the region II the function F u (t) is divided into i*fc(t)'s for t E Ik = [—kT, — {k — 1)T] 
as was done in the previous subsection. We may now study the Rondo equation, 

V 2 F k (t) = F k . x {t + T)+ v B T - (1 - k)Ax a (k > 1), (64) 

where F (t) = F : (t) and V 2 is V x with f x replaced by fc V 2 = V x - (1 - k) (X>j - 1). In 
terms of £>2, equation (j57|) becomes 

P 2 F (t) = F (t + T) + *; B T - (1 - /c)Aar n (t), (65) 

where Ax n (t) is given by (|6l|). 

6 Whcn a congested flow is realized in this model, this condition is satisfied trivially. 
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The basic technique in the last subsection may be used to solve ( |64]) with a slight 
modification. Let us define Cr&(t)'s with F k+1 (t) = F k {t) + G k +i(t + kT), which satisfy the 
equations, 

V 2 G k+1 (t) = G k (t) (k > 1), (66) 
V 2 G 1 (t)=6(e-^ t -l)=G (t) } (67) 

where 5 = (1 — k)(Axa — Axe). The boundary conditions are Gfe(O) = Gfc(O) = for 
k > 1. In solving these, we may use again the formula given in the appendix |B|. Once we 
find the series of G k {t), we obtain the function F u (t) by the relation (f44|) . Further we can 
determine the time — r by the condition (|T9|), F(—t + T) — F n (— r) + f^T = Axb- 
In the region III (t < —t), the Rondo equation take the form, 

T>iF lll (t) = F x (t + T)+ v B T + (kT 1 - l)(Ax B - Ax A ), (68) 

where F x (t + T) is F n (t + T) for -r < t + T < -r + T andF m (t + T) for t + T < -r. The 
solution must be continuously connected to F ll (t) including the first derivative at t = —r. 
Eq. (^D can be solved by the same manner as in the region II, though homogeneous 
solutions to the equation ( p8f) include the hyperbolic functions instead of the trigonometric 
functions. 

The logic at the end of section 3 may be used to get a better perspective on what we 
have discussed up to now, and it will leads us to find the a-dependence of T. First, with a 
given slope T _1 , we draw the straight line through the point S. The value of — vb and the 
coordinates of C and F are given in terms of T and the parameters in the OV-function, 

fi T a v bT . . 

Vc = T^t Vb > Axc= W^' (69) 

fiv B T + Mk' 1 - l)(Ax B - Ax A ) 

VF = i^Kf ' 

a ^T + / 1 T(/ t - 1 -l)(Ax B -Aa; A ) 

= r^If ' (70) 

and 

( f 2 T + l \ f 2 T-l 

-v B =Ui TyfT- 1 Ax A + 2T Ax B - (71) 

In solving the Rondo equation, we have introduced a parameter r determined by 
Eq. ( |I9|) . Using this r and v B expressed as Eq. (|7TD, we obtain one parameter family 
of (a-dependent) solution to the Rondo equation. Then, we find an appropriate value 
of a for a given T by the requirement that the asymptotic trajectory connects C and F: 
F m (-oo) = up. 

In Fig. [| we show trajectories from our analytic study and a computer simulation for 
the double slope model. We have chosen a particular value for a so that r = T and we 
used the Rondo function for t > — r — T to draw the part of the decelerating asymptotic 
trajectory. Clearly the Rondo function reproduces the trajectory obtained via a computer 
simulation. 
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Figure 8: Trajectories for the Double Slope Model. Thin curve shows "limit cycle" from 
simulation with a = 1.13124 when T = r = 1.58331. The analytical solution is drawn with 
thick curve. The dotted-dashed line shows the trajectory expressed by the exponential 
function with index 7 out . The dotted line shows the OV-function with f\ = 0.25, f 2 = 1.0, 
Ax A = 1.0 and Ax B = 3.0. 



By the procedure described in this subsection, we may easily obtain the remaining part 
of the asymptotic trajectory. When this is carried out, we expect that it reaches to the 
point F in the infinite past. In the following we will give another argument to support this 
expectation. The Rondo equation for t < — r — T in the region III may be solved with 
exponential functions plus a particular solution, as for the region I: the function F m (t) 
is the sum of a term linear in t and exponential functions. In the limit t — > — oo, only 
the linear term survives expressing that vehicles have the velocity for a free region, vp] the 
exponents satisfy Eq. (p9|) and have positive real parts so that exponential functions vanish 
as t — > — oo. 

It would be appropriate to explain how solutions to Eq. fl59|) are distributed on the 
complex 7-plane. Without going into the details, here we only quote features relevant to 
our arguments. There is only one solution with negative real part, it actually a real solution 
— 7i n . Other solutions have positive real parts, which are relevant when t ~ — oo. There 
is only one real solution, 7 out ; there are complex pair solutions with their real parts larger 
than 7 out . Since 7 out has the smallest positive real part, it dominates among exponential 
functions when t — > — oo. Thus for very large negative t, the function F Ill (t) may be 
approximated by the exponential functions with 7 out . In the region I, we found that the 
decelerating asymptotic trajectory is linear on the Ax-v plane owing to the exponential 
term in F (t). Similarly, the approximated F Ill (t) define a line on the Ax-v plane, starting 
the point F as shown in Fig. ||. We observe that this line is actually the tangent line to 
the trajectory at the point F. 

The absence of the flat region is directly related to our observation that the point F 
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in reached only in the infinite past. If we consider the limit to have a flat region, fi = 0, 
we only have a negative solution 7 in = —a. When we would like to consider the OV-model 
applied for a realistic situation, the relevant OV-function may be approximately realized 
by a piece-wise linear function. Since it is unlikely that the function has a flat region, the 
above feature of the double slope model must be generic. 



5 Trajectories around Cusps 

In the last section, we discussed asymptotic trajectories in models with piece-wise linear 
OV-functions. The asymptotic trajectories can be realized only when the number of vehi- 
cles becomes infinite. In computer simulations in Refs.|l], 0, a finite number of vehicles 
run around a circuit; a vehicle goes through all the free and congested regions in a finite 
time. The Rondo equation probably has solutions even for such situations. Though we 
have not worked out how to obtain entire trajectories for such vehicles, we are able to 
discuss parts of the trajectories around the points C or F. We are going to discuss this 
subject in this section. 

To make our explanation concrete, we take a trajectory around the end point C. In 
a congested region, all the vehicles have almost the same velocity. When a vehicle is 
about to reach a congested region, its velocity would be slightly different from vc and the 
function F(t) may be written as F(t) = vet + We take a linear approximation for an 
OV-function V(Ax) 7 around the point (Axc,vq), 

V{Ax) ~ f c (Ax - Ax c ) + v c . (72) 

Since {v B + v c )T = Ax c , the Rondo equation (H) becomes a linear difference-differential 
equation for £(t), 

— + &) = /o(e(* + r) -£(*))■ (73) 

a 

For the ansatz £(t) = e 7 *, we find an equation for the exponent 7, 

- + 7 = /c(e 7T -l). (74) 
a 

As long as fc < there are two real solutions: the negative, — 7; n , and the positive one, 

7out-Q 

Trajectories considered here cross (at t — 0) the OV-function at points slightly different 
from C: we denote their coordinates by (Ax,v) = (Ax c + Sv/ fc,vc + Sv). We find the 
solution for F(t) as, 

F{t) = v c t + 7out 7in 5v. (75) 

TinOoutwin Tout) 

We may find Ax n (t) and v n (t) from this solution 

a A , (a - 7in)7oute- 7in * + (a + 7out)7ine 7outt 5v 

Ax n {t) = Ax c H -, ■ : - r , (76) 

a(7in + 7out) Jc 



Here as an OV-function we have in mind a smooth, but not necessarily a piece-wise linear, function. 
8 We ignored complex solutions in this approximation since the real of those solutions are larger than 
Tout- See the discussion in the last section. 
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v n (t) = v c + 



7oute 



-lint 



+ 7ine 



7out* 



-5v. 



Tin Tout 

By eliminating the time t, we find the equation for the trajectory around the point C 
af c (Ax - Ax c ) - (a - 7in)0 - v c )V in \(a + 7out)0 - v c ) - af c (Ax - Axq) 1 



(77) 



1. (7* 



This curve shown in Fig. ^ has two asymptotes through C with slopes afc/(a — 7 in ) and 
a fc/( a + Tout); which correspond to the decelerating and accelerating asymptotic trajecto- 
ries, respectively 




Accelerating 
Asymptotic Trajectory 



Figure 9: A trajectory in linear approximation for a vehicle passing through the congested 
region with a finite size. It has two asymptotes, accelerating and decelerating asymptotic 
trajectories. The time development is indicated with arrows. 



VXAxc) 




V(Axc) 



Figure 10: The flow diagram around the point C obtained by the linear approximation. 
Two asymptotes are asymptotic trajectories. 



curves are 



Two asymptotes of (|78|) divide the Ax-v plane into four areas. In Fig. |T~0| , 
shown for solutions to Eq. ( [73| ) with various initial conditions: the point C is a saddle 
point. The linear analysis applies to the point F as well, so curves in the left-lower region 
would describe the behavior of vehicles close to a free region. 
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The condition, (v(t) —vq)/vc « 1, helps us to evaluate tc, the time interval a vehicle 
would spend around the point C; 

t 2 «— ^— (79) 

7in7out <™ 

The time tc is related to the length of a congested region L c and the number of vehicles 
in this region N c as follows, 

N c = y, L c = N c Ax c = Axjf (80) 
Therefore the size of a congested region is larger for smaller 8v. 



6 Summary and Discussion 

We have investigated the repetitive pattern formation observed in a computer simulation 
from an analytical point of view. The Rondo approach was proposed to describe the 
repetitive pattern. In addition to a and V(Ax), the difference-differential equation for the 
Rondo function F(t) contains two macroscopic parameters, T and vb, which specify the 
motion of a global pattern. In this paper we mainly paid our attention to Rondo functions 
for asymptotic trajectories. The Rondo equation was solved for three simple models with 
piece-wise linear OV-functions. In order to determine the a-dependence of T and vb, we 
gave analytic expressions for Rondo functions. We would like to emphasize that the concept 
of asymptotic trajectory plays a key role to determine the a-dependence of T and vb- As a 
first step to understand more realistic situations, we have studied some trajectories around 
cusps. 

In the following we discuss three questions related to the Rondo approach: (1) exten- 
sion to OV-models with asymmetric OV-functions; (2) more on realistic trajectories; (3) 
possibility to find the Rondo function forward in time. 

Studying more realistic models, we may encounter an asymmetric OV-function. We ex- 
plain how our procedure developed in this paper may be extended to such situations. Even 
with an asymmetric OV-function, we may define the concept of accelerating and decelerat- 
ing asymptotic trajectories. When the symmetry is absent, accelerating and decelerating 
asymptotic trajectories are not related each other and must be found independently. The 
condition that both of them share the same end points C and F will determine the ex- 
dependence of T. 

Even though we have mainly studied asymptotic trajectories, the Rondo equation itself 
must be applicable for any repetitive motion of vehicles. On the other hand, as we have 
observed for piece-wise linear models, the difference between asymptotic trajectories and 
the results of computer simulations are very small. Therefore whether we would like to 
obtain a realistic trajectory out of the Rondo approach or not is very much dependent of 
our purpose. 

The Rondo equation gives us a functional relation which may be written as, 

F(t + T) = V[F(t), F(t), F(t); a, v B , T], (81) 

with three parameters a, v B and T. When we know the function F(t) for the time interval 
t E [t ,t + T], there are two ways to use the above equation: (1) substituting this on the 
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r.h.s., we find F(t) for t G [to + T, to + 2T]; (2) the same information may be used on the 
l.h.s. to have a differential equation for F(t) on the interval t G [to — T,t }. Although the 
former sounds much easier, we have been able to use the Rondo equation only in the latter 
manner. Here we would explain the reason why it was so. 

Now let us consider, for concrete, a decelerating asymptotic trajectory in the double 
slope model. In order to use Eq. fl8l]) in the approach (1), we need the Rondo function 
describing a part of the asymptotic trajectory coming out of F; the rest of the asymptotic 
trajectory may be obtained just by differentiating the initial function repeatedly. So the 
initial Rondo function is of vital importance. 

Since the OV-function is symmetric, the transcendental equation (|59|) may be used to 
find the initial Rondo function. Let us remember how solutions are distributed. There 
are only two real solutions, — 7i n < and 7 out > 0, and infinitely many complex solutions 
whose real parts are larger than 7 out . The initial Rondo function is a linear combination of 
infinitely many exponential functions with Re (7) > 0. Therefore it contains infinitely many 
coefficients, which must be determined to give an asymptotic trajectory with properties 
described in section 3. To find an asymptotic trajectory in this way, we probably need 
some new techniques. 

In this paper we have considered the repetitive pattern in traffic flow. Such repetitive 
structure is also observed in various phenomena, and we believe that our approach may be 
helpful to understand them. 

Appendix 

A Asymptotic trajectory is symmetric 

In this paper we have used the following property of an asymptotic trajectory: (Axq,vc) 
and (Axf,vf) are at symmetric positions for OV-function symmetric around a point S. 
Here we would like to give a proof of the above claim for an OV-function symmetric with 
respect to the point S. 

Suppose an OV-function function and a sensitivity are given. Our Rondo equation 
contains two parameters T and Vb] 



If we could find a solution F(t) for the equation, this means that, for the OV-function and 
the sensitivity a, corresponding pattern with the delay T and the backward velocity v B 
may be realized. 

The OV-function is taken to be an odd function around the point S(Axs,fs)- This 
assumption is expressed with an odd function W(— x) = — W(x) as follows, 



F(t) + F(t) = V(F(t + T) — F{t) + v B T). 



(82) 



a 



V(Ax) =v s + W(Ax - Ax s ). 



(83) 



By putting this form to the Rondo equation, it now looks like, 



1 * 



F(t) + F(t) =v s + W(F(t + T) - F(t) + v B T - Ax s ). 



(84) 



a 
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We assume that a solution to the Rondo equation has been found. It gives a trajectory 
on the Ax-v plane, whose coordinate we denote as (Axi(t), Vi(t)). They are expressed with 
the function F(t) as follows, 

v 1 (t)=F(t), Axi(t) = F(t + T) -F(t) +v B T. (85) 

By using the fact that W is odd, it is easily shown that (Ax2{t) , V2(t)) given below 
satisfies the Rondo equation as well. 

v 2 (t) = F(t), Ax 2 (t) = F(t + T)-F(t)+v' B T, (86) 

where F(t) = 2v s t - F(t) and v' B T = -2v s T - v B T + 2Ax s . 

It is also easy to see that (Axi(t),Vi(t)) and (Aa^i), Vi(t)) are symmetric with respect 
to the point S. Therefore if the former defines a trajectory from a free to a congested region, 
the latter defines that for opposite direction. 

Here we emphasize that two trajectories have different backward velocities but with 
the same delay time T. 

In computer simulations, we observe that a pattern of a congested flow is characterized 
with two parameters T and v B ; both regions, connecting free to congested or congested to 
free, moves with the same backward velocity v B . So two trajectories connecting free and 
congested regions must have the same parameters. This must be also true for an asymptotic 
trajectory. Therefore v' B must be equal to v B itself. This implies that two trajectories 
discussed above form a closed trajectory. Thus we may conclude the following: 1) solutions 
expressed by F(t) and F(t) satisfy the Rondo equation with the same parameters T and 
v B ; 2) the two points on the OV-function connected by trajectories are symmetric with 
respect to the point S; 3) the straight line through the two points includes the point S. 



B General formula for Step-by-Step method 

In the following we will give a general formula for the second order linear differential 
equation, 

JHsH G * (i) = G *- i(t) 

G*(0) = 0, G k (0) = (k>l) (87) 

G? o (0) = 0, G (0) = -aS (88) 
The solution for the Eq. (|88D is 

G (t) = 6 (e- at - l) (89) 
In terms of the function gk{0) defined as follows, 

G k+X (t) = G k (t) + — ( e-^g k {ut) (k > 0) (90) 



23 



the initial conditions at t = are expressed as 

9k(0) = 0, ^(0) = 0, (k > 1) 

g (0) = 0, g' (0) = 1. (91) 
From the Eqs. (jS7|), (jSS|) and u 2 = af — a 2 /4, the equation to determine g k (9) is 



&{0) + gk{0) = g^iO) (k > 1) 

#)+9o(^) = (92) 



The initial value problem with (51) and (^) may be solved by the functions 

M-™±jM&- l)m ™-{-') (93) 



Here 

/ sin 9 \ J sin (m : even) 
I cos J 1 cos (m : odd) 



(94) 



and [r] is the maximum integer which does not exceed r. We give functions for k — 0, 1, 2, 3 
explicitly. 

0o (0) = sinfl, 

X (0) = l(sm9 -9cos9) , 



g 2 (9) = I (3sm9-39cos9-9 2 sm9 



gz [Q) = (15 sin 9- 15flcosfl-6# 2 sinfl + fl 3 cos#J . (95) 

In this article, we also use another series of solutions h k (t) of equations fl9~2"| ) with initial 
condition; 

h k (0) = 0, ^(0) = 0, (k > 1) 

/i (0) = 1, h' Q (0) = 0. (96) 
It is easily shown that h k {9) is given by 

h k {9) = g' k {9) = ^9k-x(e), (97) 

where the second equality is valid only for k > 1. We also give h k (9) for k = 0,1,2,3 
explicitly. 

ho(9) = cos 9, 

h x (9) = -9sin9, 
2 



h 2 (9) = I (9 sin 9-9 2 cos 9) , 

8 ^ ' 

h 3 (9) = — (39 sin 9 - 3# 2 cos 9 - 9 s sin 0) , (91 
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